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Abstract 

We propose a new approach to the study of the interquark potential 
in Lattice Gauge Theories. Instead of looking at the expectation value 
of Polyakov loop correlators we study the modifications induced in the 
chromoelectric flux by the presence of the Polyakov loops. In abelian 
LGTs, thanks to duality, this study can be performed in a very efficient 
way, allowing to reach high precision at a reasonable CPU cost. The 
major advantage of this strategy is that it allows to eliminate the 
dominant effective string correction to the interquark potential (the 
Liischer term) thus giving an unique opportunity to test higher order 
corrections. Performing a set of simulations in the 3d gauge Ising model 
we were thus able to precisely identify and measure both the quartic 
and the sextic effective string corrections to the interquark potential. 
While the quartic term perfectly agrees with the Nambu-Goto one the 
sextic term is definitely different. Our result seems to disagree with 
the recent proof of the universality of the sextic correction. We discuss 
a few possible explanations of this disagreement. 



1 Introduction 



An increasing amount of numerical results supports the idea that the flux 
tube joining the quark-antiquark pair in the confining regime of a gauge 
theory should be described by an effective string theory. Starting from the 
seminal paper by Liischer, Symanzik and Weisz [1] many theoretical efforts 
have been devoted in these thirty years to understand the nature of such 
string, constrain its action and predict its effects on the interquark potential. 
At the same time, the impressive development, both from the theoretical and 
from the computational point of view of Lattice Gauge Theories (LGTs) 
provided a better and better numerical laboratory to test these predictions 
(for un updated review of lattice results see: [2j). One of the most interesting 
recent results in this context is that the first few terms in the perturbative 
expansion of the string contribution to the interquark potential should be 
universal and thus can be evaluated assuming the simplest possible string 
action, i.e. the Nambu-Goto string. Universality up to the quartic order was 
proved a few years ago by Luscher and Weisz [3] and recently extended to the 
sixth order by Aharony and collaborators [IHS] (see also [7]). These results 
are based on a set of assumptions (which we shall discuss in more detail 
below) and it would be of great importance to test them with a numerical 
simulation. This is a very difficult task for at least two reasons: 

• In the standard "zero temperature" interquark potential, higher order 
corrections are proportional to higher powers in (where R is the 
interquark distance) and are thus visible only at very short distance 
where the effective string picture breaks down and spurious effects 
(due for instance to boundary terms) and perturbative contributions 
become important. 

• The dominant string correction (the so called Luscher term) may 
shadow the subleading terms 

In order to overcome these two problems we propose the following strat- 
egy. First, instead of working at zero temperature we shall study the in- 
terquark potential at finite temperature (just below the deconfinement tran- 
sition). It is well known that in this regime the string corrections are pro- 
portional to R and act as a temperature dependent renormalization of the 
string tension. In this regime higher order corrections correspond to higher 
powers of T and can be observed much better than in the zero temperature 
limit. Second, in order to eliminate the dominant Luscher term we shall 
not measure directly the interquark potential, but shall instead study the 
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changes induced in the flux configuration by the presence of the Polyakov 
loops. We shall show below that as a consequence of this choice the Liischer 
term vanishes. This makes this observable an unique tool to explore higher 
order corrections. Another reason of interest of this approach is that it is 
deeply related to another important issue of the effective string description 
of LGTs, i.e. the study of the flux tube thickness. This observable has re- 
cently attracted much interest [8j,|9j. In particular it has been shown that in 
the high temperature regime in which we are presently interested the width 
of the flux tube increases linearly with the interquark distance and not log- 
arithmicaly as one would naively expect [8l[9] . This linear increase is related 
to the linear increase in the flux energy that we observe here. In both cases 
the slope is temperature dependent and contains information on the higher 
order effective string corrections in which we are intersted. 

This paper is organized as follows. We shall introduce our approach and 
define the fiux observables in which we are interested in sect. 2 . Sect. 3 will 
be devoted to a detailed discussion of the effective string predictions for 
these observables. In sect. 4 we shall compare these predictions with a set of 
high precision simulations performed in the 3d gauge Ising model while the 
last section will be devoted to some concluding remarks and a discussion of 
the interplay of our results with the universality of the sixth order correction 
recently discussed by Aharony and collaborators in [IHS]. 

2 Flux density in presence of the Polyakov loops. 

The lattice operator which measures the flux through a plaquette p in pres- 
ence of two Polyakov loops P, P' is: 



where Up is, as usual, the trace of the orderd product of link variables along 
the plaquette p. Different possible orientations of the plaquette measure 
different components of the flux. We shall be interested in the following in 
the mean flux density, i.e. the sum of (j){p] P, P') over all the positions and 
orientations of the plaquettes, normalized to the number of plaquettes of 
the lattice. Due to translational invariance this quantity will depend only 
on the distance R between the two loops and on the inverse temperature L 
(i.e. the lattice size in the compactifled direction). Let us deflne 
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where Np denotes the number of plaquettes of the lattice. We shall denote in 
the following with Ng and L the number of lattice spacings in the spacelike 
and timelike directions respectively so that in d = 2 + 1 we have Np = 3N^L 
It is easy to see that if we define: 

Z{R,L) = {PHR)Pm = i j;pt(^)p(o) e'^E.c/, (3) 

c 

where the partition function in eq.Q is given as usual by 

Z = ^e^Sp^- (4) 

c 

and J2c denotes the sum over all the configurations, then {^{R,L)) can be 
written as: 

mR,L)) = ^^^logZ{R,L) . (5) 

If we neglect for the moment effective string corrections and keep only the 
area term in Z, i.e. Z{L, R) ~ e~'^^^ we find a linearly rising behaviour for 

{<^{R,L)) = aR (6) 

with an angular coefficient 

L da 

'' = -N-pTp 

which does not depend on the finite temperature 1/L (which cancels out in 
the ratio L/Np) 

If we now look at the effective string corrections to this result we imme- 
diately observe a very interesting feature of {^{R,L)) i.e. the fact that the 
dominant term (the so called Liischer term) vanishes. 

In fact one of the main assumptions of the effective string approach to 
LGTs is that the only dimensional parameter which appears in the actiorE] 
is the string tension a. Therefore observables like Z{R,L) may depend on 
/3 only through the string tension a. This means that in the effective string 
framework we can convert the /3 derivative into a derivative with respect of 
a. 

However it is well known [1] that the dominant correction must be scale 
invariant and thus a independent. This means that in eq.Q this term 
disappears and the effective string correction to {^{R,L)) starts from the 
first subleading term: the "quartic" correction that we shall discuss below. 



^In principle one could also add to the action a "boundary term" proportional to a 
dimensional parameter b. This term vanishes at the first order in (2+1) dimensions [3] 
but it may appear at higher orders. We shall discuss below this issue. 
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3 Effective string corrections 

In this section we shall first briefly recall a few general results on the Nambu- 
Goto string and on its perturbative expansion (a detailed derivation can be 
found in [3] and [TO]) and then use these results to discuss, in the framework 
of the Nambu-Goto effective string, how the flux distribution is modified by 
the presence of the Polyakov loops. 

3.1 The Nambu-Goto string 

The action of the Nambu-Goto model is proportional to the area of the 
string world-sheet: 



where g is the determinant of the two-dimensional metric induced on the 
world-sheet by the embedding in and a is the string tension, which 
appears as the only parameter of the effective model. 

Eq. ([8]) is invariant with respect to reparametrization and Weyl trans- 
formations, and a possible choice for quantization of the effective model is 
the "physical gauge" in which g becomes a function of the transverse dis- 
placements of the string world-sheet only. The latter (which can be denoted 
as X*(r, <f)) are required to satisfy the boundary conditions relevant to the 
problem. In the present case, the effective string world-sheet associated with 
a two-point Polyakov loop correlation function obeys periodic b.c. in the 
compactified direction and Dirichlet b.c. along the interquark axis direction: 



This gauge fixing implicitly assumes that the world-sheet surface is a single- 
valued function of (r, ?), i.e. overhangs, cuts, or disconnected parts are 
excluded. It is well known that rotational symmetry of this model is broken 
at the quantum level because of the Weyl anomaly, unless d = 26. However 
this anomaly is known to vanish at large distances |llj . and this suggests 
to use the "physical gauge" for an IR, effective string description also for 



Here and in the following, we restrict our attention to the d = 2 + 1 
case, which is particularly simple, as there is only one transverse degree of 
freedom (X). In the physical gauge, eq. ([8]) takes the form: 




X*(r + L,0=X*(r,0; 



X\t,0) = X'{t,R) =0 



(9) 



d / 26. 




(10) 
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This action describes a non-renormalizable, interacting QFT in two di- 
mensions; the associated partition function is expected to encode an effective 
description for this sector of the gauge theory, providing a prediction for the 
VEV of the two-point Polyakov loop correlation function: 

(Pt(i?)p(0)) = Z{L,R) = J PXe-^[^] . (11) 

3.2 Perturbative expansion of the partition function 

The effective model can be studied perturbatively, expanding the square 
root appearing in eq. (|10p in powers of the dimensionless quantity {aRL)~^: 
this approach is expected to be consistent with the fact that the effective 
theory holds in the IR regime of the confined phase. We find: 

-aRL ry ( , , 



Z(L,i;)=.--.Z,.^l + -^ + ^ + ...J (12) 

where the indices in F4 and recaU the fact that they are obtained from 
the quartic and sextic terms in the expansion of the Nambu-Goto action 
respectively. 

The leading order of this expansion: Zi corresponds to the partition 
function of a free boson in two dimensions: 

(13) 



7]{iu) 



where u = and r] denotes the Dedekind function: 



r]{r)=q^^llil-qn ; g = e^-- , (14) 

n=l 

This correction is universal and does not depend on the string tension a. 
Thus, as anticipated in the introduction, it will not contribute to {^{R,L)). 

The next to leading corrections was evaluated in |12j and can be written 
as a combination of Eisenstein functions. 

where E2 and are defined as: 



E2{t) = l-24^ai(n)(z" (16) 



n=l 
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E^{t) = l + 240j^a3(n)g" (17) 

n=l 

q = e^-- , (18) 



where o"i(n) and (73(71) are, respectively, the sum of all divisors of n (includ- 
ing 1 and n), and the sum of their cubes. 

This expression looks complicated, but it simplifies in the large R/L 
limit (i.e. g — )• 1) in which we are interested. By performing a modular 
transformation we can rewrite the above expression in terms of g = e"^'^^ . 
In the large R/L limit we may neglect all higher order terms in the q ex- 
pansiorll and after some algebra we end up with the following contribution 
to {<^{R,L)): 

We see that as a consequence of the flux tube fluctuations a 1/R term appear 
in the mean flux density and that the linerly rising term of eq.® gets a 
correction proportional to 1/L^ and thus becomes temperature dependent. 

Higher order terms, like Fq, could be derived, in principle, along the 
same lines but, if one is interested in the Nambu-Goto action only, they 
can be more simply obtained by expanding in power of l/{aLR) the exact 
solution obtained a few years ago by Liischer and Weisz [3j and then derived 
in the covariant formalism in [13j. 

In d = 2 + 1 dimensions one finds that Z(R, L) is given by a tower of Kq 
Bessel functions [3t[T3]: 

00 

Z{R, L) = (P(0, 0)P(0, R)) = CnKoiEnR). (20) 

n=0 

where En are the closed string energy levels: 



Svr 

En = aLU + 



and a denotes the zero temperature string tension. 



1/2 

(211 



^We shall make this statement more rigorous in the section devoted to the numerical 
simulations. 
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It is easy to see that in the large R hmit only the lowest state (n = 0) 
survive^ and we end up with a single Kq function: 



with 



and 



lim {PiO,0)P{0,R)) =coKo{EoR). (22) 

R— s>oo 



/ vr \i/2 



C^lf- (24, 

From eas.(|22l23land [2i|) we find 

dlogZ{R,L) _( I ^ j,KdEo\ da 

dl3 ~ [la^'^Ko da J dp ^ 

where we used the fact that Z{R, L) is a function of j3 only through the 
string tension a 

Using eg. ([23]) we find 

dE^ ^ L{l-x/2) 

da ^/T^ ' 
where we shall denote from now on 



Recalling that 



and using the asymptotic expansion of the modified Bessel functions we 
obtain 

^<iiEoR) ^ . , ^ L_ + (2R) 

Kq{EoR) ^ 2EoR 8iEoR)^ ^ ^ ^ 

Collecting all the terms together we find 

mR,L)) = a(RA{x) + Bix) + ^ + ..) (29) 



^This statement is equivalent to the choice of neglecting terms proportional to q that 
we made in the previous section. 
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where we defined: 

AW = S^ (30) 

and a was defined in eq.© which we report here for the particular case of 
a (2+1) dimensional lattice for which Np = 3N^L: 

While a depends on the details of the gauge theory and in particular 
on the /3 dependence of the string tension a, the functions A{x),B{x) and 
C{x) ... encode the information on the effective string model. The particular 
form listed above for these functions corresponds to the Nambu-Goto model, 
however according to |3,-'i] the first two orders in the perturbative expansion 
in power of l/o" (i.e., with our conventions, the terms proportional to l/c^ 
and 1/c^) powers of x) should be universal. In order to identify these terms, 
let us expand in powers of x these functions. 

Expanding in x we find 

l + y + y + -J (34) 
B{x) = ^(^ + ^ + ..) (35) 



Cix) = ^rZ7^{l + x + ...) (36) 



o-L V 4 4 
1 

(^)2 

which at the first order coincide as expected with the corrections ob- 
tained with the previous approach. In addition we find the next to leading 
corrections of order 



3.3 Boundary corrections 

The analysis of [3] excluded the presence of a boundary term at the first 
order (see eq. (2.4) of [3J) but could not exclude the presence of boundary 
terms at higher orders. These higher order boundary terms were recently 
studied in great detail in [6]. Looking at their result it is easy to see that 
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they cannot contribute to the term proportional to R in the expectation 
value of {^{R,L)), but they may affect the following terms (constant and 
proprotional to 1/R)- For this reason in the following we shall concentrate 
only on the linear correction. 

4 Numerical simulations 

4.1 Details on the computational setting 

We compared our predictions with a set of high precision simulations in 
the 3d gauge Ising model. We performed our simulations using the same 
methods discussed in Oni]. We used duality to map the Polyakov loops 
correlator into the partition function of a 3d Ising spin model in which 
we changed the sign of the coupling of all the links dual to the surface 
bordered by the two Polyakov loops. We then estimated L)) by simply 

evaluating the mean energy in presence of these frustrated links. We chose 
periodic boundary conditions in the original gauge Ising model. These b.c. 
are mapped by duality into a mixture of periodic and antiperiodic b.c. in the 
dual spin model. However we always chose A'^, large enough to eliminated 
any contribution from the antiperiodic sectors (which are depressed by terms 
proportional to g^'^^'^ or e^'^^*). It is clear that duality plays a crucial 
role in this derivation and for this reason we consider the approach discussed 
in this paper as particularly suited for abelian gauge theories, even if in 
principle, given enough computational power, it could be extended also to 
non-abelian models. 

Since, as discussed above, we shall be interested in the following only 
in the term proportional to R in {^{R,L)), (which is the only one which is 
not affected by boundary corrections) we shall neglect in the following the 
disconnected component (Up) of {^{R,L)). 

4.2 Evaluation of a for the 3d gauge Ising model 
4.2.1 Naive analysis 

Assuming the scaling behaviour 

a(/3)=^7,(/3,-/3)2- (37) 

we find 

_ Ida ^ 2u 
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which would lead to a very simple ed elegant form for a 



However this is not enough for our purposes. In order to estimate higher 
order effective string corrections we need to evaluate the flux density with 
an uncertainty smaller than 1% and thus it is mandatory to have a rigorous 
control of all the possible sources of systematic error up to this threshold. 
We have in particular four potential sources of systematic errors: the deter- 
mination of the critical point /3c and of the critical index i/, the determination 
of the zero temperature string tension a and the presence of scaling correc- 
tions which could affect the derivative Let us addres these point in more 
detail 



4.2.2 Scaling corrections and critical parameters. 

The most accurate result for the critical temperature of the 3d spin Ising 
model is (3c,s = 0.22165455(5) [15] from which we can immediately obtain /3c 
using the duality relations. The most accurate values given in the literature 
for the critical index u are u = 0.63012(16) [16j from the analysis of high 
temperature series expansions and v = 0.63002(10) [l7j from a recent finite 
size scaling analysis of Monte Carlo data. We shall adopt this last value in 
the following analysis. For both /3c and v the available precision is much 
smaller than our statistical errors thus we may neglect this possible source 
of systematic error. While there are no precise enough results on the scaling 
behaviour of the string tension we can extract the information that we need 
from a very precise study appeared a few years ago [18] on the scaling 
behaviour of the dual observable: the interface tension of the 3d Ising spin 
model. Let us define as /3s the dual of the gauge coupling /3 and accordingly 
as t = /3s — /3c,s the reduced temperature of the 3d spin model. 

In the neighbourhood of the transition, the interface tension behaves as: 

cr{(3)=crct^'' x{l + at^ + bt + ct^^ + df + ...) , (40) 

The most accurate result for the subleading exponent is = 0.832(6) [17] . 
from which 6 = vuj = 0.5241(33). For 9' we may use the result 0' = 
1.05(7) [H]. 

Since 20 « 0' ~ 1, we may assume a simplified ansatz for the scaling 
behaviour: 

(7(/3) = Uct^ X (1 + at^ + ht) . (41) 
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Using a set of extensive simulations of the interface tension in [18] the 
authors were able to extract rather precise estimates for the free parameters 
of the above ansatz: Uc = 10.083(8), a = -0.479(26) and h = -2.12(19) 

We are interested in the the value of: 



Let us rewrite it as: 



1 da 



usmg 



Ida 
adp 



1 dPs da 1 djig da 
a d/3 dPs a df3 dt 



1 



d/3 sinh2/3 



we find, 



1 da 



1 



a dp sinh2/3 V t 



2u 



1 H f + -^^ —t 

2u 2u 



(42) 

(43) 
(44) 

(45) 



where we made the approximation 20 ~ 1 which allows us to combine 
the correction proportional to with the one proportional to b. 

In the range of values of /3 that we studied this correction turns out to 
be larger than our statistical errors and thus must be taken into account 

Using the above scaling function we may also obtain very precise esti- 
mates for the string tension for the values of /3 that we studied 



4.3 Results and comparison with the effective string predic- 
tions 

We per fomed simulations for three values of /3: /3 = 0.75180, 0.756427, 0.758266 
which were chosen because for these values the deconfinement temperature 
is known with very high precision and coincides with l/T^ = = 8, 12, 16 
respectively [20]. In particular, we studied several values of L in the /3 = 
0.75180 case and a single value of L for /3 = 0.756427, 0.758266 to test the 
scaling behaviour of the corrections. For each value of L we studied a wide 
set of values of the interquark distance in the range < R < . For each 
value of R and L we used 10^ iterations to thermalize the lattice and then 
performed 10^ measures. Errors were evaluated with the standard blocking 
and jacknife techniques. 

Let us see in detail the analysis of the /3 = 0.75180 data. 
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4.3.1 Results for (3 = 0.75180 

The results of the simulations are reported in tab{T]For this value of /3, the 
value of a obtained from the scaling function eq. (j4ip is cr = 0.0105255(11) 
(see |18]). It is interesting to notice that in this case we also have an estimate 
for sigma directly obtained in the gauge theory: a = 0.0105241(15) from [21] . 
The fact that the two values agree within the errors represents a strong 
crosscheck of the procedure that we adopted to extract a. In the following 
we shall use the first estimate which is slightly more precise. Using this 
value and the scaling analysis discussed in the previous section we obtain 
for a the following prediction a = 2.792 10~^. To test this prediction and 
at the same time the form of the effective string correction we carefully 
selected two sets of values of L. In the first set we chose three values of L. 
L = 16, 20, 24 for which the variable x (defined in eq. ()27p ) was tuned so as 
to make negligible within the errors the sextic term 2:^/8, but not negligible 
the quartic one x^/8 in eq. (l34p . In the second set we chose three values of L: 
L = 10, 11, 12 for which instead the sextic term was definitely not negligible 
within the errors. 





L = 10 


L = 11 


L = 12 


L = 16 


L = 20 


L = 24 


R = 8 


N/D 


0.913241(5) 


0.912778(4) 


0.912136(3) 


0.912022(3) 


0.911996(3) 


R = 10 


N/D 


0.913300(5) 


N/D 


N/D 


N/D 


N/D 


R = 12 


0.914175(7) 


0.913357(5) 


0.912896(4) 


0.912251(3) 


N/D 


N/D 


R = 14 


N/D 


0.913421(5) 


N/D 


N/D 


N/D 


N/D 


R = 16 


0.914303(7) 


0.913486(5) 


0.913013(4) 


0.912362(3) 


0.912249(3) 


0.912225(3) 


R = 20 


0.914442(7) 


0.913619(5) 


0.913138(4) 


0.912484(3) 


N/D 


N/D 


R = 24 


0.914563(7) 


0.913733(5) 


0.913256(4) 


0.912593(3) 


0.912472(3) 


0.912450(3) 


R = 32 


0.914790(8) 


0.913969(5) 


0.913488(4) 


0.912823(3) 


0.912697(3) 


0.912664(3) 


R = 40 


0.915052(8) 


0.914194(5) 


0.913730(4) 


0.913052(3) 


0.912921(3) 


0.912893(3) 


R = 48 


0.915279(8) 


0.914442(6) 


0.913958(5) 


0.913284(3) 


0.913144(3) 


0.913121(4) 


R = 56 


0.915524(9) 


0.914688(6) 


0.914188(5) 


0.913506(3) 


0.913371(3) 


0.913344(4) 


R = 64 


0.915765(9) 


0.914906(6) 


0.914423(5) 


0.913734(3) 


0.913601(3) 


0.913568(4) 



Table 1: Values of $(i?, L) for p = 0.75180 



For each value of L we fitted the data for $(i2, L) according to the law 

<^{R,L) = a{L)R + b{L) + c{L)/R (46) 

where the term c{L)/R was introduced only for the second set of values of 
L because it was always compatible with zero in the first set. The results of 
the fits are reported in tabl2J 
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a{L) 


c(L) 




L=10 


0.913840(9) 


3.016(21) e-05 


-0.000213(55) 


1.29 


L=ll 


0.913033(10) 


2.944(20) e-05 


-0.000257(77) 


1.87 


L=12 


0.912560(7) 


2.917(13) e-05 


-0.000145(68) 


0.84 


L=16 


0.911914(3) 


2.8472(62) e-05 


N/D 


0.81 


L=20 


0.911795(5) 


2.8163(86) e-05 


N/D 


0.66 


L=24 


0.911773(3) 


2.8025(78) e-05 


N/D 


1.37 



Table 2: Results of the fits to eq.® for /? = 0.75180 



In the following we shall concentrate on the values of a{L) which as 
discussed above is the only one which is not affected by boundary corrections. 

The values of a{L) extracted from the fits are reported in tablH and 
plotted in fig.l. We analyzed these data in two steps. First we fitted the 
first three values of a{L) (those corresponding to L = 16, 20, 24) with the 
law: 

a(L) = a(l +7x2) ^^j^ 

we found a = 2.7918(17) 10"^ and 7 = 0.132(7) with a very good x^- Both 
these values nicely agree with the predictions a = 2.792 10~^ and 7 = 
1/8. Notice that if we would neglect the scaling correction discussed in the 
previous section and use eq. (j39|) to estimate a we would find a = 2.807 10~^ 
in complete disagreement with the numerical result. This gives an idea of 
the precision of the our estimates and of the need of carefully controlling 
scaling corrections. 

We then fitted the whole set of data (i.e. including also L = 10, 11, 12) 
with the law: 

a{L) = a(l + 7^2 + 6x^) (48) 

we found again a good reduced x^- Xr — 0-^^ the following best fit 
results: 

a = 2.796(5) 10~^ 7 = 0.127(25), 5 = -0.051(27) 

The first two values agree again very well with the predictions but the coeffi- 
cient of the sextic correction, which should he 6 = 1/8 completely disagrees. 
This can be well appreciated looking at fig.l where we plotted the quartic 
correction (dashed curve), the sextic correction according to the Nambu- 
Goto effective action (dashed-dotted curve) and that corresponding to the 
best fit value of the parameter 6 (dotted curve). 
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Figure 1: Plot of a A(L) as a function of ^ = j-. The continuous line is the pre- 
diction for a (i.e. without effective string corrections. The other curves correspond, 
from top to bottom to the whole Namhu-Goto prediction, ea ll30\) . the truncation at 
the sextic order, the truncation at the quartic order and finally the dashed dotted 
line corresponds to the best fit choice (5 = — 0.05j of the sextic coefficient (see the 
text). The points are the results of the simulations in the 3d gauge Ising model 
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L 


a{L) 








10 


3.017(21) 


2.792 


3.137 


3.481 


11 


2.943(20) 


2.792 


3.027 


3.222 


12 


2.917(13) 


2.792 


2.959 


3.074 


16 


2.847(9) 


2.792 


2.845 


2.865 


20 


2.816(9) 


2.792 


2.814 


2.819 


24 


2.802(8) 


2.792 


2.802 


2.804 



Table 3: Values of the coefficient a{L) for various values of L. In the second 
column we list the results of the simulations performed at /3 = 0.75180. In the 
following columns we report the prediction for a{L) = aA{x) at the zeroth, first and 
second order in the expansion in x 

4.3.2 Results for /? = 0.756427 and /3 = 0.758266 

In order to test the scaling behaviour of the effective string corrections we 
also simulated a value of L corresponding to T/Tc = 0.8 (i.e. the highest 
value of T among those discussed in the previous section) for two other values 
of j3 nearer to the critical point. More precisely we performed simulations 
at /3 = 0.756427 on a lattice of size 15 x 192^ for values of R in the range 
12 < E < 84 and at /? = 0.758266 on a lattice of size 20 x 256^ for values 
of R in the range IQ < R < 128. In order to keep under control the 
statistical uncertainty we perfomed in these cases 15 x 10^ and 5 x 10^ 
measures respectively. We report the results of the simulatioins in tabS] 
and [5j Moreover, in order to eliminate spurious effects due to a possible 
breaking of the string picture for small values of R we decided to fit only 
the data in the range R > ^ . With this choice the 1/R term in the fit was 
always compatible with zero. Thus we decided to eliminate it from the fit 
and perform only two parameters fits (constant plus linear term). 

The results are reported in tabE] and figl2] where we have plotted the 
function A{x) defined in eq. (j3l|) as a function of the finite temperature T/Tc. 
Since both quantities are adimensional we may plot different values of f3 on 
the same graph. In figiH we also reported the Nambu-Goto prediction for 
A{x) (dotted line) together with the value truncated at the fourth order 
(dashed line) and at the sixth order (dashed dotted line). 

It is easy to see that all the points lie below the fourth order line (i.e. 
they are compatible with a negative value of 5) and show a slight tendency 
to move toward the value (5 = as /3 moves toward the continuum limit. In 
all the cases the data are incompatible with 5 = 1/8. This test seems to 
suggest that the deviations that we observe from the universal sixth order 
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L = 15 


R = 12 


0.928413(4) 


R = 18 


0.928496(5) 


R = 24 


0.928572(5) 


R = 30 


0.928638(5) 


R = 36 


0.928716(5) 


R = 42 


0.928785(5) 


R =48 


0.928841(5) 


R = 54 


0.928922(5) 


R = 60 


0.928976(5) 


R = 66 


0.929051(5) 


R = 72 


0.929120(6) 


R = 78 


0.929187(6) 


R =84 


0.929257(6) 


R = 90 


0.929335(7) 


R = 96 


0.929387(6) 



Table 4: Values of L) for (3 = 0.756427 

term is not due to the effect of an irrelevant operator but, at least in the 
range of values that we studied, seems to hold up to the continuum limit. 

5 Conclusions 

In this paper we discussed a new strategy to test effective string predictions. 
Instead of looking at the interquark potential we studied the flux distribu- 
tion in presence of the quark- antiquark pair. The major consequence of this 
choice is that the leading effective string correction (the Liischer term) dis- 
appears thus allowing to disentangle higher order corrections. We apply the 
method to the 3d gauge Ising model and find that while the first subleading 
correction (the one related to quartic terms in the string Lagrangian) agrees 
very well with the one predicted by the Nambo-Goto model, the subsequent 
term (due to sixth order terms in the Lagrangian) is definitely different and 
seems to have the opposite sign or at most to be compatible with zero. This 
observation agrees with a set of similar results obtained in these last years 
in the 3d gauge Ising model looking at various different physical observables 
ranging from the interquark potential to the effective string width to the 
interface free energy of the dual spin Ising model [9l[10l[18l[21] . These tests 
were all supporting, even if only at a qualitative level, a value for the sextic 
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0.2 0.4 0.6 

T/Tc 



0.8 



Figure 2: Plot of A{L) as a function of As in fig.l the continuous line is 
the prediction for a (i.e. without effective string corrections. The other curves 
correspond, from top to bottom to the whole Nambu-Goto prediction, ea \30\) . the 
truncation at the sextic order and the truncation at the quartic order. The points 
are the results of the simulations (see text). Stars denote the data at (3 — 0.75180, 
the black square is the value for j3 — 0.756427, T = 0.8Tc and and the open square 
the point at P ^ 0.758266 
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L = 20 


R = 16 


0.934812(8) 


R = 24 


0.934874(8) 


R = 32 


0.934907(8) 


R = 40 


0.934962(7) 


R = 48 


0.935019(8) 


R = 56 


0.935045(8) 


R = 64 


0.935090(6) 


R = 72 


0.935149(8) 


R = 80 


0.935189(8) 


R = 88 


0.935240(9) 


R = 96 


0.935279(8) 


R = 104 


0.935344(6) 


R = 112 


0.935398(11) 


R = 120 


0.935429(12) 


R = 128 


0.935477(10) 



Table 5: Values of^{R,L) for (3 = 0.758266 



correction different from the Nambu-Goto one and fully compatible with the 
one that we find here. 

Our result is rather puzzling in view of the recent proof of the universality 
of the Nambu-Goto action up to this order in d = 3 [IHS]. We see a few 
possible explanations for this disagreement. 

• It could be due to an irrelevant operator and thus it would disappear 
in the continuum limit. While we cannot exlude in principle this pos- 
sibility the scaling test discussed at the end of the previous section 
suggests that this is not the case. 

• The effect that we observe could be due to a compensation between 
the sixth order and all the remaining higher order terms which are not 
constrained by universality. While this is certainly possible, it would 
require a rather strong fine tuning of the coefficients of all the higher 
order terms. If this is the explanation, then it would be important 
to understand the reason of such a cancellation, which most probably 
doesn't occur by chance but would be the signature of some underlying 
constraint or symmetry. 

• It could be due to the fact that 3d gauge Ising model does not admit 
a weakly coupled effective string description. This limit is a basic as- 
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/? 


T/Te 


A{x) 








0.751800 


0.333 


1.0036(29) 


1.0037 


1.0044 


1.0045 


0.751800 


0.4 


1.0086(32) 


1.0077 


1.0097 


1.0102 


0.751800 


0.5 


1.020(3) 


1.0189 


1.0262 


1.0304 


0.751800 


0.667 


1.052(4) 


1.0597 


1.1009 


1.1773 


0.751800 


0.727 


1.067(6) 


1.0845 


1.1540 


1.3967 


0.751800 


0.8 


1.090(8) 


1.1237 


1.2468 


7.0472 


0.756427 


0.8 


1.082(7) 


1.1237 


1.2468 


7.0472 


0.758266 


0.8 


1.118(13) 


1.1237 


1.2468 


7.0472 



Table 6: Values of the coefficient A{x) for various values of (3 and T jT^. In the 
third column we list the results of the simulations. In the following two columns 
we report the prediction for A{x) truncated at the first and second order in the 
expansion in x. In the last column the all orders result for A(x). 

sumption of [3llll[13] and amounts to ask that the partition function 
of the string describing a particular surface (say the cylindric surface 
connecting the two Polyakov loops) can be written as a sum of single 
string states propagating along the surface. The fact that this argu- 
ment could be relevant for the 3d gauge Ising model is also supported 
by the intuitive observation that a gauge theory based on the Z2 group 
(and even more the percolation gauge theory discussed in |22) ) is in- 
deed the farthest possible choice with respect to the large limit of 
the SU(N) gauge theory which is known to behave as a weakly coupled 
string theory. If this is the reason it would be interesting to under- 
stand why the effective string description works instead so well up to 
the quartic order. 
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